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[^Abstract: We prove that for spacetimes solving the Einstein-Maxwell (EM) equations, 
the electromagnetic field contributes at highest order to the nonlinear memory effect of grav- 
itational waves. In [3] D. Christodoulou showed that gravitational waves have a nonlinear 
memory. He discussed how this effect can be measured as a permanent displacement of test 
masses in a laser interferometer gravitational wave detector. Christodoulou derived a precise 
formula for this permanent displacement in the Einstein vacuum (EV) case. We prove in 
Theorem [6] that for the EM equations this permanent displacement exhibits a term coming 
from the electromagnetic field. This term is at the same highest order as the purely gravita- 
tional term that governs the EV situation. On the other hand, in Chapter [3j we show that 
to leading order, the presence of the electromagnetic field does not change the instantaneous 
displacement of the test masses. Following the method introduced by D. Christodoulou in [5] 
and asymptotics derived by N. Zipser in [8] and [9], we investigate gravitational radiation at 
null infinity in spacetimes solving the EM equations. We study the Bondi mass loss formula 
at null infinity derived in [9]. We show that the mass loss formula from [9j is compatible 
with the one in Bondi coordinates obtained in [3]. And we observe that the presence of 
the electromagnetic field increases the total energy radiated to infinity up to leading order. 
Moreover, we compute the limit of the area radius at null infinity in Theorem [7j 
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1 Introduction and Main Results 

In this paper we investigate the null asymptotics for spacetimes solving the Einstein-Maxwell 
(EM) equations, compute the radiated energy, derive limits at null infinity and compare them 
with the Einstein vacuum (EV) case. We show that the presence of the electromagnetic field 
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does not affect the instantaneous displacement of the test masses of a laser interferometer de- 
tector at leading order, as it only comes in at lower order. But the electromagnetic field does 
contribute to the nonlinear effect of the displacement of the test masses. The EM case gives 
us a wonderful opportunity to observe mass loss and also to measure gravitational radiation. 
It is crucial to understand fully the behavior of the gravitational field also when other fields 
are present and to investigate their interplay. The only way to achieve this, is to compute 
the null asymptotics of the spacetimes. 

A major goal of mathematical General Relativity (GR) and astrophysics is to precisely de- 
scribe and finally observe gravitational radiation, one of the predictions of GR. In order to 
do so, one has to study the null asymptotical limits of the spacetimes for typical sources. 
Among the latter we find binary neutron stars and binary black hole mergers. In these pro- 
cesses typically mass and momenta are radiated away in form of gravitational waves. Bondi, 
van der Burg and Metzner studied these in [3]. D. Christodoulou showed in his paper [5] that 
every gravitational-wave burst has a nonlinear memory. The insights of this work are based 
on the precise description of null infinity obtained by D. Christodoulou and S. Klainerman 
in [7J. Among the many pioneering results they derived the Bondi mass loss formula. This 
is all in the regime of the Einstein vacuum equations. Then N. Zipser studied the Einstein- 
Maxwell equations in [SJ, [9] and computed limits along the lines of [7J for this case. She 
derived a Bondi mass loss formula, where in addition to the one obtained by Christodoulou 
and Klainerman, a component of the electromagnetic field comes in. Thus the mass radiated 
away goes into the gravitational and the electromagnetic field. Here, we rely on the methods 
introduced in [7J, used in [8], [9] and by one of the present authors in [I], [2]. There is a 
large literature about gravitational radiation. However, in the present paper, we only give 
the references which are relevant to our investigations. 

The main results of this paper are the following. We first recall the Bondi mass loss for- 
mula obtained in [9] for spacetimes solving the EM equations. 

d -M (u) = ±- / ( \E\ 2 + I \A F \ 2 ) dfic 



du y ' 8vr J S 2 V 2 1 1 y 7 

Compared to the formula obtained in [7J for spacetimes solving the EV equations, we have an 
additional term, |Af| 2 , from the electromagnetic field. Furthermore, we compare the above 
mass loss formula to the corresponding formula in Bondi coordinates [4j 

^M(ui) =-J (W) 2 + (d w d) 2 + \{X 2 + Y 2 )^ 

and show that the two formulae agree. 

As shown in the work of Christodoulou [5 J , S + — X~ is the term which governs the permanent 
displacement of test particles. Using this fact, Christodoulou shows that the gravitational 
field has a non-linear "memory" which can be detected by a gravitational-wave experiment 
in a spacetime solving the EV equations. We will describe this experiment in the last section 



as well. In section 2.3 of our paper, we study the permanent displacement formula for un- 
charged test particles of the same gravitational-wave experiment in a spacetime solving the 
EM equations and show that the electromagnetic field contributes to the nonlinear effect. 
We first obtain Theorem [6] which determines S + — £~ in the EM case. From Theorem |6j 
we observe that the electromagnetic field changes the leading order term of the permanent 
displacement of test particles. Then in the last section, we study in details a gravitational 
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wave experiment for our findings. We observe that the electromagnetic field does not enter 
the leading order term of the Jacobi equation. As a result, to leading order, it does not 
change the instantaneous displacement of test particles. But the electromagnetic field does 
contribute at highest order to the nonlinear effect of the permanent displacement of test 
masses. Furthermore, in Theorem [JJ we compute the limit of the area radius r on any null 
hypersurface C u as t goes to oo and show that the result coincides with the one obtained in 
[7] for EV. 

We follow the method introduced by Christodoulou in [5] to study the effect of gravitational 
waves. The treatment is based on the asymptotic behavior of the gravitational field obtained 
at null and spatial infinity. These rigorous asymptotics allow us to study the structure of the 
spacetimes at null infinity. The spacetime is foliated by a time function t and by an optical 
function u. The corresponding lapse functions are denoted by eft and a. Each level set of u, 
C u , is an outgoing null hypersurface and each level set of t, Fit is a maximal spacelike hyper- 
surface. We pick a suitable pair of normal vectors along the null hypersurface. The flow along 
these vector fields generates a family of diffeomorphisms (j) u of S 2 . We use cj) u to pull back 
tensor fields in our spacetime. This allows us to study their limit at null infinity along the 
null hypersurface C u . Then we study the effect of gravitational waves by taking the limit as u 
goes to ±00. Christodoulou in [5] gives a complete explanation of the structure at null infinity. 

The methods introduced in [7] , used in [8] , [9] , reveal the structure of the null asymptotics of 
our spacetimes. In these works, to prove the stability result, the data is assumed to be small. 
However, as far as the study of the null asymptotics is concerned, the data can be large. We 
give a brief outline in the last part of this introduction of the methods of [TJ . 

In General Relativity the fundamental equations are the Einstein equations linking the cur- 
vature of the spacetime to its matter content. 

G>i, := Rfj, v — - g^u R = 8ir T^ v , (1) 

(setting G = c = 1), fi,u = 0,1,2,3. G^ u denotes the Einstein tensor, R^ u is the Ricci 
curvature tensor, R the scalar curvature tensor, g the metric tensor and is the energy- 
momentum tensor. 



Here, we are discussing the Einstein-Maxwell equations. This means that T„„ on the right 
hand side of ([TJ is the stress-energy tensor of the electromagnetic field. The twice contracted 
Bianchi identities imply that 

D V G^ = . (2) 

This is equivalent to the following equation, namely, that the divergence of the stress-energy 
tensor of the electromagnetic field vanishes: 

= (3) 

with 

Tpis = ^ziFffFvp ~ ~^9puF pa F p ) (4) 

where F denotes the electromagnetic field. Note that F is an antisymmetric covariant 2- 
tensor. As T^ v is trace-free, the Einstein equations ([TJ take the form 

Rfiu = . (5) 
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We find that the scalar curvature is identically zero. We write the Einstein-Maxwell (EM) 
equations as 





= 


D a F a p 


= 


D a *F afj 


= 0. 



(6) 
(7) 
(8) 

As a consequence of the Maxwell equations, we have 

OF = (9) 
where □ is the de Rham Laplacian with respect to the metric g. 

We split the Riemannian curvature Ra^s into its traceless part, namely the Weyl tensor 
Wapiti, and a part including the spacetime Ricci curvature R a p and spacetime scalar curva- 
ture R: 

Ra^S = Wafi^S + -(g a -yR/38 + 9(5sRa-y ~ Q^RaS ~ QaSR^) 

-g idajgps - g a sg^,) R • (10) 

One can then define the Bel-Robinson tensor 

Qa^s = w apia w p p s a + *w ana *W// . (11) 

The Bianchi equations for the Weyl tensor in the presence of an electromagnetic field read 

D a W a ^ 5 = l(D y Rp S - DsRfr) . (12) 

In [7j Christodoulou and Klainerman derived the asymptotic behavior in the case of strongly 
asymptotically flat initial data of the following type. 

Definition 1 We define a strongly asymptotically flat initial data set in the sense of ^ 
(studied by Christodoulou and Klainerman) to be an initial data set (H,g,k), where g and 
k are sufficiently smooth and there exists a coordinate system (x 1 ,^ 2 ,^ 3 ) defined in a neigh- 
bourhood of infinity such that, as r = (Ylf=i( xl ) 2 )^ ~^ °°; 9ij an d hj ar ^ : 

2M . „ 3 „ , 

9ij = (1 + — ) kj + o 4 (r"5) (13) 

hj = 03 (r~i) , (14) 

with M denoting the mass. 

Under a smallness condition on the initial data, Christodoulou and Klainerman proved in [7] 
that this can be extended uniquely to a smooth, globally hyperbolic and geodesically com- 
plete spacetime solving the EV equations. The resulting spacetime is globally asymptotically 
flat. Together with the existence and uniqueness theorem comes a precise description of the 
asymptotic behavior of the spacetime. While the smallness condition was imposed in order to 
ensure completeness, the results about the behavior at null infinity are largely independent 
of the smallness. The decay behavior of the components of the Weyl tensor are given below. 
And the limits at null infinity of the relevant quantities are given in section [2j 
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In [7J as well as in [8], [9] and [T], [2] the Weyl tensor W in (M,g) is decomposed with 
respect to the null frame e4,e3,e2,ei. That is, e^ and e3 form a null pair which is supple- 
mented by e A , A = 1, 2, a local frame field for St jU = H t n C u . Given this null pair, e3 and 
e4, we can define the tensor of projection from the tangent space of M to that of St,u- 

We decompose the second fundamental form fcy of Ht into 

fcjVJV = S (15) 

= £yl (16) 

kAB = r\AB (17) 

where TV" is the unit normal vector of St >u in Ht. Let T be the future-directed unit normal to 
H t . We define 

0ab = (V A N,e B ). 

The Ricci coefficients of the null standard frame T — N,T+N,e2,e± are given by the following 



Xab = Gab - Vab (18) 

X AB = ~^ab - vab (19) 

& = ^iM-a-^a (20) 

= r 1 ^-^ (21) 

G = r'^ + eA (22) 

1/ = -<r 1 f N </>+s (23) 

1/ = ^f^ + S (24) 

J = ^-a^y^a (25) 
We use x> X) etc f° r the Ricci coefficients of the null frame a _1 (T — N), a(T + N), e2, ei. 

Definition 2 We define the null components of W as follows: 

»,u(W) = n/^^ej^ (26) 

^(W0 = InfWwseZelei (27) 

= ^ 1*W 4 el e 5 4 (28) 

a (W) = i *W a(97<5 ef e* (29) 

/3m W = ^ n/ el el e 5 4 (30) 

(W) = n/ n/ W P7<t5 e] e| . (31) 



The estimates in [7J yield the decay behavior: 

a(W) = O (r- 1 rl 5 ) 

p(W) = O (r~ 3 ) 

a(W) = O (r- 3 tZ^) 

a(W), f3{W) = o (r~i) 



5 



where t 2 = 1 + u 2 and r(t,u) is the area radius of the surface St }U - 

In [8], [9], Zipser works with the same conditions on the metric, second fundamental form and 
curvature, in addition she imposes a decay condition on the electromagnetic field F, namely 

F\ H = 03 (r~a) • (32) 
The null components of the electromagnetic field are written as 

(33) 



Fas = a{F) A F AA = a (F) A 



F 34 = 2p (F) Ftz = <J{F). 

The corresponding null decomposition {a (*F) , a (*F) , p (*F) , a (*i 7 )} of *F is given by 

a(*F) A = -a(F) B e BA a (* F) A = a {F) B e BA 
p(*F) = a(F) a(*F) = -p(F) (34) 

where the Hodge dual of a tensor u tangent to St )U , is defined by 

* B 
u A = e A u B . 

The estimates in [8], [9] yield the decay behavior: 

a(F) = O {r~ l rlh 

p(F),a(F) = 0(r- 2 rZ h ) 
a(F) = o (r~^). 

One of the main difficulties in [7j is that a general spacetime has no symmetries and thus 
does not have suitable vectorfields to construct integral conserved quantities. To overcome 
this difficulty, Christodoulou and Klainerman use the 'closeness' of their spacetimes to the 
Minkowski spacetime and construct quasi-conformal vector fields. The main step is carried 
out within a bootstrap argument in the 'last slice', namely in a spacelike hypersurface which 
is a level set of the time function t. First, the authors foliate the spacetime by functions 
t and u near the initial slice. From the foliations, one constructs vectorfields that are al- 
most Killing. Combining these vectorfields with the Bel-Robinson tensor, one obtains local 
estimates for the Weyl curvature tensor W and the electromagnetic field F. With these 
estimates, one constructs a new optical function which is defined on a larger domain in the 
spacetime. Then, following a continuity argument, one obtains a smooth, globally hyperbolic 
and geodesically complete spacetime solving the Einstein equations. The resulting spacetime 
is globally asymptotically flat, satisfying the above decay properties. 

Acknowledgment. We thank Demetrios Christodoulou for helpful discussions, critcal read- 
ing of our paper and constructive feedback. 

2 Null Asymptotics 

2.1 Asymptotic Behavior and Bondi Mass 

We need precise data at null infinity. In particular, we have to know the Bondi mass and the 
asymptotic behavior of the components of the curvature and the electromagnetic field. Zipser 
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described them in [8], [9] following the discussion in Chapter 17 of [7], making changes as 
necessary due to the presence of the electromagnetic field. The parameters of the foliations 
and the components of the Weyl tensor behave exactly as in [7|. That is, the following holds: 
along the null hypersurfaces C u as t — > oo, it is 



and 



lim <f> = 1, lim a = 1 (35) 

Cy,,t—>-oo C u ,t->oo 



lim {rtrx) = 2, lim (rtrx) = —2 (36) 

C u ,t— >oo C u ,t— >oo 



Furthermore, we let 



H = lim 

O u ,t— >oo 



(r\tr X ' - 2 r ) \ . (37) 



From the existence theorem of jS], [9], Zipser makes the following conclusions, which are 
generalizations of conclusions 17.0.1 through 17.0.4 in [7J. 

Following the convention in [7] and [9J, the pointwise norms | | of the tensors on S 2 re- 
late to the metric 7, which is the limit of the induced metrics on St, u rescaled by r~ 2 for each 
u as t — > 00 . 

Theorem 1 On any null hypersurface C u , the normalized curvature components ra(W), 
r 2 /3(W), r 5 p(W), r 3 a(W), ra(F), r 2 p(F), r 2 a (F) have limits as t — )• 00, in particular 



lim ra (W) = A w (u, ■) , lim r 2 (3 (W) = B w («, •) 

C u ,t— >oo C u ,t— >oo 

lim r 3 p(W) = iV («,-)> lim r 3 a(W) = Q w (u,-) 

C u ,t— >oo C u ,t— >oo 



lim ra(F) = Ap (u, •) , 

Cujt— > OO 

lim rVt-F) = P F (u,-), lim r 2 cr (F) = Q F (u, •) 

C u ,t— >oo C u ,t— >oo 

wit/i j4vf a symmetric traceless covariant 2-tensor, Bw and Ap 1- forms and Pw, Qw> Pf> 
Qf functions on S 2 depending on u. The following decay properties hold: 

\AwM\ < C(l + \u\)- 5/2 \B w (u,-)\<C(l + \u\r 3/2 

\P w (u,-)-P w (u)\ < (1 + lnl)- 1 / 2 \Q w (u,-)-Q w {u)\ < (1 + H)- 1 / 2 
\A F (u,-)\ < C(l + M)- 3/2 

\PfM\ < (1 + lnl)" 1 / 2 \QfM\ < (1 + ln))- 1 / 2 



and 



lim Pvv (u) = 0, lim Q w (u) = 0. 

M— >— OO 11— > — OO 



The existence of the limits in the conclusion follows from the estimates in the existence 
theorem of [9] (i.e. 0). 

Theorem 2 On the null hypersurface C u , the normalized shear r 2 x' has limit as t — > 00: 

lim r 2 x=S(it,-) 

U u> t— >oo 

with S feeing a symmetric traceless covariant 2-tensor on S 2 depending on u. 
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The proof is the same as in [7] because the propagation equation stays unaltered 

—— = -trxxAB ~ a{W )ab- 



Theorem 3 On any null hypersurface C u , the limit of rrj exists as t — > oo, that is 



lim rrj = S (it, 
C U} t— >oo 



with 5 fremg a symmetric traceless 2-covariant tensor on S 2 depending on u and having the 
decay property 

|s(«,.)|. <c(i + Hr 3 / 2 . 



Further, it is 



lim r# = lim r%' = S 

C u ,t— >oo 2 C u ,t— 5>oo — 



as weZZ as 

OS _,_ 

oe _ _i 

5u ~ ~4 



(38) 

A w . (39) 



Zipser proves this result as conclusion 3 in [9]. The argument is along the lines of the proof 
of conclusion 17.0.3 in [7J. 

Zipser follows [7J to derive the Bondi mass formula by calculating a propagation equation for 
the Hawking mass enclosed by a 2-surface St 7 u- The Hawking mass is defined as 

m(t,u)= T - + / tr X trx). (40) 

\ J St.u J 

Let 

M = + • X " P 0*0 - ^ {p 2 (F) + a 2 (F)) . (41) 

With respect to the /-pair, one has the null structure equations 

—= + - trx tr K = _ 2 ^+2|C| 2 

dJ ^ + \^xf = -\x\ 2 ~\a{F)\\ 

One computes, 

d - trxtrx + trx {trxtrx) = - 2jJ,trx + 2tr x | C | 



ds 



2 

2 . I /-T71\|2 



-trx\x\ -trx\a(F)\' 



thus 

d_ 

dt 



I trxtrx = -2 a<j)fxtrx 

J St,u J St,u 

+ / a<p(-trx\x\ 2 -trx\a(F)\ 2 +2trx\([ 



(42) 
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Using the Gauss equation 

K = -\tr X trx + \x ■ X ~ P (W) - \ (p 2 (F) + a 2 (F)) , 

one derives 

jx = -di/vQ + K + -trxtrx- (43) 



By the Gauss-Bonnet formula and formulas ( |41[ ), (43), conclude that 

/ /£ = / fe-X-pW-kp 2 (^) + cx 2 (F))) 

J St u J St.u V / 



(44) 



Moreover, by 



47r ( 1 + T7T I tr Xtrx | = — m. 



<i r- 



-r = -Satrx, 
dt 2 Y A ' 



and (42), (44), it is 



9 r 



m(t,u) = -— / (afarx - <f>atrx) (J< (45) 



St,u 



r P 1 1 1 1 

J St.u V 

Note that i<C + \tr X trx = 0(r~ 3 ), /i = 0(r -3 ). From the asymptotic behavior of the 
right-hand side of (45), it follows 

— m (t,u) = O (r~ 2 ) . 

This means that m(t,u) has a limit for any fixed u as i — >• oo, namely the Bondi mass of 
the null hypersurface C u . As in [7J, it is denoted by M (it). The terms appearing due to the 
presence of the electromagnetic field are shown to decay fast enough so that the mass decays 
at the same rate as in [7J. In particular, 

m(t,u) = M{u) + (r -1 ) 

as t — > oo on C u . 

Following [7J, Zipser calculates a Bondi mass loss formula by considering 

-m(t,u) 

with 

Q i r r 

—m (t, u) = -atrOm + — — / a N N (i + tr9/j) . 

As I = a' 1 (T + N) and l = a(T — N), 

1 



a(V Nl± + tr6^) = -a (D^fi + trxp) 



2 
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and 



D 4 ^ + trx/i = 0(r" 4 ) 

T> 3 fi + trxjJ- = ~^trx\x\ 2 - -trx\a(F)\ 2 + 0(r~ 4 ) 



Thus, 



d r 

Tu m{t ' u) = ^ 



j s ^x{\x\ 2 + \\a(F)\ 2 ^ +0(r- 1 ). 



Following [7J , Zipser uses the following facts to derive the Bondi mass loss formula: the metric 
7 = 4**11 ( r_2 7) converges to the standard metric 7 of the unit sphere S 2 as t — > 00 for each 
u (44 u ^ s a diffeomorphism from S 2 to St lU ), moreover \trx converges to 1, and r% converges 
to — 2S. This yields 



<9it w 8vr 7 5 2 V 2 1 1 / 'V 

The right-hand side of this expression is positive and integrable in it. Thus, M (it) is a non- 
decreasing function of it and has finite limits M (—00) for it — > —00 and M (00) for it — > 00. 
Further, Zipser concludes from (44) that M (—00) = 0, and M (00) is the total mass. 

Theorem 4 The Hawking mass m(t,u) tends to the Bondi mass M (u) as t — >■ 00 on any 
null hyper surf ace C u . That is, 

m(t,u) = M(u) + 0(r _1 ). 
and M (v) verifies the Bondi mass loss formula 

~M (u) = — [ (\E\ 2 + 1 I A F \ 2 ) dfio 
8u K ' Sir J S 2 V ' 2 1 1 J ^ 

with d/j,o being the area element of the standard unit sphere S 2 . 

We see that in the Bondi mass loss formula the limiting term Ap of the electromagnetic field 
comes in. At this point, let us compare this with the Bondi mass loss formula obtained in 
|7J (p. 499): J^M(u) = J s2 \E\ 2 d/j,a. In fact, the electromagnetic field contributes to the 

change of the Bondi mass by J s2 |^4f| 2 dfio. 

The decay behavior of Ap is the same as for S. See theorems [T] and [3j Similarly as in 
[7], [5] for the Einstein vacuum case, we can define now the new function 

Then is the total energy radiated to infinity in a given direction per unit solid angle. Thus 
the integrand in ( 46 ) is proportional to the power radiated to infinity at a given retarded time 
u, in a given direction, per unit area on S 2 (per unit solid angle). Already in [5] Christodoulou 
tells us how to adapt the formula for F when matter radiation is present, that is also in the 
EM case. 



In the next two subsections, we also need the following theorem for H . 
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Theorem 5 The function H satisfies 



H = (48) 



Proof: In the EV case, equation (47) is proved in conclusion 17.0.5 of [7J where one uses the 
fact that 

^Ntrx' + \x' = 0(r~ 3 ). 

In the EM case, it is easy to see that the additional terms involving the electromagnetic field 
are also 0(r -3 ). Thus equation (47) is still true in the EM case. 



In the EV case, equation ( |48| ) is proved in lemma 17.0.1 in [7J. In the proof, we need to 
show that r 2 5 converges to 2M(u). From Proposition 4.4.4 in [7], we have 

4irr 3 5 = j du'(f arO ■ i) k(5 — 5) — ra~ l ^f7 a ■ e + r{divk)^) 



o 

- «o J St,u * 
Following the proof of lemma 17.0.1 in [7], we see that 

/ ar6-rj--K(5-5)-ra- 1 fa-e = r [ \E\ 2 dfio + O(l) 
Js t ,u 2 Js 2 7 

Moreover, in the EM case, due to the constraint equation, we have 

(divk) N = R 0N = 2F p F Np . (49) 



Using the equation ( 49 ) , we see that 

/ r(divk) N = - \A F \ 2 dfio + O(l) 
JSt.u 2 Js 2 7 

since the leading term of Fqa and Fna are both \o.{F)a- As a result, we can still conclude 
that 

r8 = —p: I r — m(t,u) + 0{r~ l ) 
r 2 Ju du 

The rest of the proof follows easily from the proof of lemma 17.0.1 in [7J. 
2.2 Compare result with mass loss in Bondi coordinates 

In this subsection, we compare the mass loss formula obtained in [8], [9], cited above in The- 
orem [4j with the mass loss formula in Bondi coordinates. Bondi coordinates are first defined 
by Bondi, van der Burg and Metzner in [3] for axially symmetric vacuum spacetimes. The 
main motivation for such coordinates is to study gravitational radiation at null infinity. The 
form of the metric is chosen such that many computations are simplified at null infinity. As 
a result, one can derive many useful theorems and formulae assuming the existence of such 
coordinates. In particular, the Bondi mass loss formula in Bondi coordinates is first derived 
for axially symmetric vacuum spacetimes in [3] and is generalized to the Einstein-Maxwell 
case in [3]. However, for a given spacetime, it is hard to tell whether such coordinates exist. 
On the other hand, spacetimes studied in j8] are obtained by evolving small initial data on a 
spacelike hypersurface by EM equations. In particular, it is not clear whether all spacetimes 
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studied in [5J admit such coordinates. Here we show that for the leading term in the mass 
loss formula, the two different coordinate systems give the same result. 



The mass loss formula from [Hj, [9] is: 

^-M (u) = — I f |~| 2 + - \A F A duo. 
du y 1 8vr J S 2 V 2 1 1 / 7 

Let us recall the mass loss formula and asymptotic expansion for solutions to the Einstein- 
Maxwell equations in Bondi coordinates [5]. The line segment is: 

-UVdw 2 - 2Udwdr + o ab (dx a + W a dw)(dx b + W b dw) a, b = 2, 3 

with the electromagnetic field given by a skew-symmetric two tensor Fu U . We have the 
following asymptotics for the line segment. 

V = 1 - — + 0(r~ 2 ), U = 1 + 0(r- 2 ) and W a = 0(r" 2 ). 
r 

_ ( r 2 + 2cr + • • • — 2dr sin 6 + • • • 

<?ab ~ y _ 2dr sind + ... sin 2 Q( r 2 _ 2cr ) + . . 

We also need the following asymptotics for 

F we = X + 0{r~ l ) and F W(k = Y sin# + 0{r~ l ). 

as well as 

F ra = 0(r- 2 ), F ab = 0(1) and F wr = 0( r - 2 ). 

Level sets of w, C w , are outgoing null hypersurfaces. Each C w is then foliated by level sets 
of r, S w>r . Hence, it is natural to consider the following pair of null vectors normal to S wr 

dw dx c 2 dr dr 

since 64 is a natural choice of null vector on C w and 

lim (e 3 ,e 4 ) = -1. 

r— >oo 

Let M(w) be the Bondi mass. It is given by 

M(w) = — / mdu,o. 
V ' 8ttJ S 2 ^ 



The mass loss formula reads 



j KS„c) 2 + (fl„<i) 2 + \(X* + Y 2 )) dp. 



To show that the two mass loss formulae agree, we prove that 

\E\ 2 = (d w c) 2 + {d w d) 2 and \A F \ 2 = X 2 + Y 2 . 
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First we compute 



// d 9 _ d ^ c d V d d 

— dx a ' dx b dw dx c 2 dr ' ai^ dx b 

^ dx c 2dr ,ab dr^ ab dw ab dx c 

d _„„ 9 V „„ 3 „„„ 5 , 



/_ — yy c — — — — r r ,— + r u- 

K dw dx c 2 dr' ab dr ab dw 



The Christoffel symbols are 

Kb= -\g wr d r g ah = \d r g ah + 0{l) 

T ab = ^9 wr {db9wa + d a g w b - d w g ab ) + ^g rr {-d r g a b) + ^g rc (dbg ca + d a g cb - d c g ab ) 
= \d w g a b - \d r g a b + 0(1) 



As a result, one finds 

One can easily see that up to O(l) terms, d w a a b is traceless and <9 r <7 a 6 has zero traceless part. 
As a result, 

|H| 2 = (d w c) 2 + (<9^) 2 . 

For the second equality, we use the expression for e% and the asymptotics for F^ v . A direct 
computation shows that 

l^l 2 = X 2 + Y 2 . 
2.3 Permanent Displacement Formula 

The permanent displacement of the test masses of a laser interferometer gravitational-wave 
detector is governed by £ + — £~. Christodoulou showed in [5] how this works. We discuss the 
corresponding wave experiment in the EM case in section [3j In the following, we are going to 
state and prove a theorem for £ + — £~ in the EM case. We point out that the final formula 
- even though in its form identical to the one obtained by Christodoulou and Klainerman in 
|7J - differs from the EV case by a contribution from the electromagnetic field. The form of 
the formula is not altered due to the fact that the corresponding extra electromagnetic terms 
cancel. However, the limiting term Ap enters the new formula nonlinearly. 



Theorem 6 Let £+(■) = lim 

u^oo £(^> •) and £ (•) — lim ll _ ! ._ 00 £(u, •). Let 

/OO 1 
(\E(u,.)\ 2 +-\A F (u,-)\ 2 )du . (50) 
-oo ^ 

Moreover, let $ be the solution with <& = on S 2 of the equation 

$=F-F . 

Then £ + — £~ is given by the following equation on S 2 : 

di/v (£+ - £-) =f $ . (51) 



13 



u — > — oo. Also, it is 



Proof: Equation (38) in theorem [3] yields that £ tends to limits S + as u — > oo and S as 

/« 
E(u)du' 
-oo 



and 



/oo 
E («')&*' . 
-oo 



When taking the limits for the Hodge system ((|57|), (|58|)) on C n as t — > oo, we will compute 
the corresponding limits for and involving . From Zipser's work [9], chapter 9.1, (9.13) 
and lemma 9, we know 

A4> = r | fj | 2 - r - | a(F) | 2 (52) 

2 —1 

r z a 1 



A*' = -ra" 1 A( | fj | 2 -| r/ | 2 ) + (a# 4 | a(F) | 2 -a|) 4 | a(F) | 2 ) (53) 

whereas in the work of Christodoulou and Klainerman [7], chapter 11.2, (11.2.2b) and 
(11.2.7b) it is 

Al' = r | fj | 2 (54) 

A*' = -ra~ 1 \( | fj | 2 -| f] | 2 ) . (55) 

We compute the following limits in the new case. 

lim ^ = * lim ^' = 

lim rV^W = fi(u,-) lim rV N ^' = n'(u, •)• (56) 

C u ,t— >oo C' u ,t— >oo 

We proceed by investigating the Hodge system for e. The Hodge system for e reads, see also 
[9], chapter 9: 

3 

dijve = -V N 5 - -tr65 + f] ■ 6 

^(a-^-e + J |a(F)| 2 -J |a(F)| 2 (57) 

qiWe = <t(W) + 0At/ . (58) 

We observe that the qiW equation coincides with the one obtained by Christodoulou and 
Klainerman in [7J, whereas the di/v equation contains the extra terms | a(F) | 2 and | a(F) | 2 
from the electromagnetic field. (See [TJ, chapter 17, ((17.0.12a), (17.0.12b)).) According to 
[9], chapter 9, one has: 

V N 5-e-fj+ - | a(F) | 2 = -2r~ 3 (V N r)p + r~ 2 V N p-r~ 2 (V N r)V N y + r~ 1 V 2 N ^ 
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-X-r]- r 1 ^-r 2 (rtr9 + a 1 A)Vat* 

-r^a^fa ■f^ + r~ 2 V N p - 2r~ 3 a~ 1 \p (59) 



with 



p = rV nQ + q + ^>' 

This differs from [7J, chapter 17, (17.0.12c) by the extra curvature term from the electromag- 
netic field. Zipser derived in [9], chapter 9, Lemma 9, 

Aq = r{pL-JL)+I (60) 
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where 



I 2 -3 I «00 



-A* 



— . fj — K fi — 2ra l y a ■ e + - \ a(F) \ 2 
and is the mass aspect function given by 

p = -p(W) -X-fj ■ 

Recall the radial decomposition of A to be = A — trOV^ — A - a^fa -f. Now, we 
obtain from the last equations that 



Aq 



V%q + tr6V N q + Aq + a'^a- fq 

—r(p — p) — r\ ■ f) — kS — 2raT x y a ■ e + - | a(F) | 2 



(61) 



We proceed as follows: Substituting first for V^vP from (61) in (59) and then the resulting 
terms from (59) in (57) yields 

di/ue = p- p + X-V~X-V + r~ l - r~ 2 V N ^>' - r~ 3 a _1 A*' + l.o.t. (62) 
qiWe = a(W) + 6Af) (63) 



Let 



E = lim (re) 

C'u,t— >O0 



We multiply equations (62) and (63) by r and take the limits on C u as t — > oo. This yields: 



ci^rl E 
di/u E 



q + s a: 



p-p + s-s-s-s 
+ a * - - n' . 



(64) 



(65) 



Then we investigate the limits as u — > +oo and u — > — oo. Considering the last equations 
for e, respectively E, and using theorems [3] and [TJ one finds that E tends to a limit E + as 
u — > +oo and to E~ as u — > — oo. 



By conclusions along the lines of [7], chapter 17, we obtain 

ci/lH (E + -E~) = 

o 

In order to compute diju (E + — E~ ) , we have to consider especially the corresponding limits 
for the terms involving and that is also Q' . 

Much like Christodoulou and Klainerman computed the formulas in lemma 17.0.2, on page 
504 of [7], we derive the new results in which the electromagnetic field term a(F), respectively 
its limit Ap, is present. To do that, we use the fact that 

P A a{F) A = --tr X a(F) A + l.o.t. (66) 
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The derivation of (66) can be found in Zipser's work [9 J on p. 351 formula (4.15). Now, 
considering ( 53 ) and using ( 66 ) , we find that 

Pa I ! 2 = ~ tr X I &{F) | 2 +l.o.t. (67) 

Using (67), (54) and (55), we deduce formulas for 'J/', £1, £1' by computing the limits (56). 
We give the formulas: 



n 



+ 0O 



EfjvW) , 1 f \A F | 2 (u',u/ 



—da; + 

22 Air J -oo ^Js 2 (l-wo/)2 - 



s 2 (l-uu')i 



j — doj 



1 



22 47TJ-00 L JS 2 



1 



2247T^-oo l is2 (l-ww') 5 



(1 — CJ0/)2 



/5 2 



(1 — OJUJ 1 ) 2 



52 (1 



- WW 2 



1 /■+ 



2i47r 



| 2 (u',u') 



00 <• Js 2 



(1 — CJU)') 2 



(1 — WW')2 

{s<?n(u -«')(( I s I 2 " I 3 I 2 («')) + j( I ^ | 2 («>') - \ A F f («' 



}• 



Straightforward calculation shows that when evaluating the difference of the limits as 

o 

u —> +00 and u —> —00 in (65), the contribution of If\ ty' and fi' comes only from terms 
in We find that Q! tends to limits r2 /+ (-) and Q!~ (•) as t — )• 00 and i — > — 00, respectively. 
Thus, we conclude 



+00 



1 , 1, 



= jf ( |H(u,-)| 2 -|H( U ,-)| 2 + -|^K-)| 2 -2l^(n,-)| 2 )^ • (68) 
Finally, we obtain 

di/u(E + -E~) = + (69) 

+00 1 1 

( - I S(tv) I 2 +1 S(«,.) I 2 - - I A F (u,-) I 2 +-| A F (u r ) p )d« . 

-00 z z 

Proceeding along the lines of [7], chapter 17, it is 

(E+-E~) = f $ (70) 
with $ being the solution of vanishing mean of 

£ $ = + jy- on s 2 . 



Accordingly, also by conclusions along the lines of [7J, chapter 17, we derive (72). To see this, 
we consider the normalized null Codazzi equation 

(di/ux)A ~ l? A tr X + e B XAB - \e A tr X = ~P{W) A - p{F)a{F) A - e AB a{F)a{F) B (71) 
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Multiply equation (71) by r 3 and take the limit as t — > oo on C u . We obtain 



di/v S = f H + E 



as in [7] p. 510, conclusion 17.0.8 since the extra terms from the electromagnetic field in (71) 

(72) 



decay fast enough. Due to equation (47), we conclude 

di/v (S + - 5T 



E + -E~ 



Thus, the theorem is proven. 



2.4 Limit for r as t — > oo on Null Hypersurface C u 

We shall use the fact that the constraint on the spacelike scalar curvature, which is given by 

R=\k\ 2 + R 00 , 

differs from the constraint in the vacuum case only by the term Roo, which is a quadratic in F. 

Building on the results of D. Christodoulou and S. Klainerman in [7] as well as the results of 
N. Zipser in [9] (i.e. [8]), we can now prove the following results. 

Theorem 7 As t — > oo we obtain on any null hypersurface C u 

r = i-2M(oo)logt + 0(l) . 
Proof: We recall from p. 503, with <fi' = <fi — 1, 

dr r- 



dt 



4>trx' 



^(i + 0')(- + (^x / --)) 

2 r r 

i + W + o( r - 2 ) 



In the last equality, we use equation ( 48 ) . 
From [9j, p. 465 



R 



1 



oo 



(\a(F)\ 2 + \a(F)\ 2 ) + p(F) 2 + a(Ff 



(73) 



with a(F),a(F), p(F),a(F) the components of the electromagnetic field. Moreover, the lapse 
equation in our situation is given by 



k \ 2 +R t 



oo J 



(74) 



We integrate the lapse equation (|74|) on H t in the interior of S t)U ' to obtain 

du' I a<P(\k\ 2 +R 00 ) . 



j Vjv0'= f du' j 

J St, u Jug J S t , 



In view of ( 73 ) and the fact that all the terms on the right hand side of ( 73 ) except a(F) are 



of lower order, we estimate 

St,u 



du 



o0(| k | 2 +- | a{F) \ 2 ) + l.o.t. 
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We see that 

/ a<t>{\ k | 2 +\ | a(F) | 2 ) -> / | H | 2 +^ | A F 
Consider the Bondi mass loss formula in theorem [4j Then, as i — > oo we conclude 

Vtv^' - 8ttM(u) = 0(r- 1 ) (75) 



' St,u 

on each C u . In view of (/)' we compute: 
u 1 /" 1 



div{r-*<t>'N) 



1 f 00 1 



/ oVat^' + (atr^ - atr6)<f/) 



4ttJ u {r{t,u')f 

i f°° i r 



where B denotes the exterior of St tU - Therefore, from (75) it follows on C u as t — > oo, 

<f/(t,u) = -2 , . ,,. „ M{u')dv! + 0(r- 2 ) = — M (oo) + 0( r - 2 ) . 
Ju {r(t,u')y r 

Thus, we obtain on any cone C u for t — > oo, 

dr 1 

— = 1 - -M(oo) + 0(r- 2 ) . (76) 
at r 

Thus, the statement of our theorem follows, which closes the proof. 



3 Wave Experiments 

We are now going to show how the results above relate to experiment. In (Sj Christodoulou 
established his breaking result on the nonlinear memory effect. The idea of the gravitational 
wave experiment and setup is given in [5 J , discussing a laser interferometer gravitational- wave 
detector. There Christodoulou explained how the theoretical result on S + — X - leads to an 
effect measurable by such detectors. This effect manifests itself in a permanent displacement 
of the test masses of the detector after a wave train has passed. In the present EM case, 
we find a result on the displacement of test masses which is twofold. Considering the Jacobi 



equation (see (86)), the highest order term remains unchanged. However, there is an extra 
term at highest order from the electromagnetic field in the formula for S + — S - , the perma- 
nent displacement of test masses, as we have shown in the proof of theorem [6] In the present 
chapter, we shall show how the electromagnetic field enters the experiment and we will derive 
results for this case. 

We will follow the lines of argumentation by Christodoulou in [5] and [6]. 
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Let us briefly review the setup of a laser interferometer experiment: Three test masses are 
suspended by pendulums of equal length. Denote by mo the reference mass, which is also 
the location of the beam splitter. For timelike scales much shorter than the period of the 
pendulums, the motion of the masses in the horizontal plane can be considered free. By laser 
interferometry the distance of the masses mi and mi from the reference mass mo is measured. 
Whenever the light travel times between mo and mi and m2, respectively, differ, this shows 
in a difference of phase of the laser light at mo- 

The masses mo, mi, m,2 move along geodesies 70, 71, 72 in spacetime. Let T be the unit 
future-directed tangent vectorfield of 70 and t the arc length along 70. For each t denote by 
H t the spacelike, geodesic hyperplane through 70 (t) orthogonal to T. 

Take (E\, E2, E3) to be an orthonormal frame for H$ at 70 (0), and parallelly propagate 
it along 70 to obtain the orthonormal frame field (T, E\, E2, E%) along 70. It follows that 
(Ei, Ei, E%) at each t is an orthonormal frame for Ht at 7o(i)- Then one assigns to a point p 
in spacetime, lying in a neighbourhood of 70, the cylindrical normal coordinates (t , X , X , X ^ . 
based on 70, if p £ H t and p = expX with X = £\ x % Ei 6 T l0 ^H t . Denote by d the distance 
of p from the center 7o(i) on H t , that is, d =| X \= -^/^(x*) 2 . The difference between the 
metric and the Minkowski metric r]^ u in these coordinates is: 

9^ - Vk* = O (R d 2 ) . (77) 

As usual, we put c = 1. Now, let r be the time scale in which the curvature varies significantly. 
Then, the displacements of the masses from their initial positions will be 0(Rt 2 ). Assume 
that 

d 

- « 1 . (78) 
r 

Then we can read off from the differences in phase of the laser light differences in distance of 



mi and m-2 from mo- Further, in view of (78) one can replace the geodesic equation for 71 



and 72 by the Jacobi equation (geodesic deviation from 70) 

d 2 x k 



dt 2 



RkTIT x l (79) 



with RkTIT = R {Ek, T, Ei,T). One is free to assume that the source is in the ^-direction. 
This was derived by Christodoulou for the EV case in [S], and can also be found in his [B]. 



We now investigate the formula ( 79 ) for the Einstein-Maxwell situation. If we assume the test 



masses not to be charged, then formula (79) stays the same, but through the EM equations 



and in view of (10) the electromagnetic field comes in. We shall see that it enters at lower 



order though. From ( 10 ) one can write 



Rkom — Wkoio + - (gkiRoo + 9oqRh — doiRw — 9koRoi) ■ (80) 
As there is from the EM equations: 

Roo = 87tToo , 

and in particular, we have 

^00 = i(| a(F) | 2 + I a(F) | 2 ) + p(F) 2 + a(F) 2 (81) 
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we can investigate the components of the Ricci curvature on the right hand side of ( 80 ) . The 
component i?oo includes the term | a(F) \ 2 . Recall that a(F) is the part of the electromag- 
netic field with worst decay behavior. However it enters as a quadratic the formula for Rqq. 



To proceed, we consider L 
vature are 



T — E3, L = T + E3. The leading components of the cur- 



QUb(W) 
<*ab(W) 



R (E A , L, E B , L) 



A AB (W) 



+ o (r 



-2\ 



(82) 
(83) 



And the leading components of the electromagnetic field are 



a A (F) 
a A (F) 



F{E A ,L) 
Aa(F) 



+ o (r 



(84) 
(85) 



In the following the kth Cartesian coordinate of the mass m A for A = 1,2 will be denoted 

k 

(A)' 

d 2 x k 



by x k , A y Then the Jacobi equation becomes 



(A) 



dt 2 



- r 1 A AB x l 



1 



(B) 



A F I 2 x l 



(B) 



+ O {r~ z ) 



that is 



d 2 x 3 



(C) 



dt 2 



2 „A 



d l x 



(C) 



dt 2 



1 _1 . B 1 

-r A AB x {D) - - 



r- 2 I Ap | 2 x\ D) + O (r- 2 ) 



(86) 

(87) 
(88) 



From the Jacobi equation (86) we see that the electromagnetic field enters on the right 
hand side at order (r~ 2 ) only. Thus, we have shown that the electromagnetic field does 
not contribute at leading order to the deviation measured by the Jacobi equation. There- 
fore, at leading order, we can rely on the results for the Einstein vacuum case, derived by 



Christodoulou in [5]. Instead of (86) he obtained 



d 2 x> 



(A) 



dt 2 



- - r 1 A AB x\ B) + 0(r 



,-2\ 



(89) 



As in [5] one obtains that in the vertical direction there is no acceleration to leading order 
(r _1 ). Initially mi and ttt-2 are at rest at equal distance do and at right angles from mo- 
This implies the following initial conditions, as t — > —00: x' 



d 5% , x B (A) 



(A) 







x 3 



(A) 







r B 
J (A) 



x"r A \ = 0. The right hand side being very small, one can substitute the initial values 
on the right hand side. Then the motion is confined to the horizontal plane. One has to 
leading order: 



..A 

X (B) 



1 



do A AB 



One obtains 



(B) 



(*) 



A AB (it) du 



(90) 



(91) 
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In the following, let us revisit the result (95) from [5]. In view of equation (39), i.e. 



dE 
du 



and Hmi^i^^ S = we obtain 



/' 

j — i 



A A b (u) du 



and 



As 



x: 



(B) 



do 



Iab (t) 



A 



(92) 



(93) 



for u — > oo, the test masses return to rest after the passage of the gravitational 



wave. Taking into account (38), i.e. ^ = — S, and integrating again: 



x: 



(B) 



(0 



,d . 



i^ AB (t) - S- 



(94) 



The limit t — > oo is taken and it follows that the test masses experience permanent displace- 
ments. Thus S + — £~ is equivalent to an overall displacement of the test masses: 



A x J 



,do, 
r 



(B) ~ - '< ■.. K^AB ~ S Ab) • ( 95 ) 

The right hand side of (|95|) includes terms from the electromagnetic field at highest order 



as given in our theorem 



6l Even though the form of (95) is as in the EV case investigated 



by Christodoulou in [H] and [B], the nonlinear contribution from the electromagnetic field is 
present in — 
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